Parity Doubling and the S Parameter Below the Conformal Window 
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We describe a lattice simulation of the masses and decay constants of the lowest-lying vector and axial 
resonances, and the electroweak S parameter, in an SU(3) gauge theory with Nf = 2 and 6 fermions 
in the fundamental representation. The spectrum becomes more parity doubled and the S parameter per 
electroweak doublet decreases when Nf is increased from 2 to 6, motivating study of these trends as Nf 
is increased further, toward the critical value for transition from confinement to infrared conformality. 
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Introduction In a recent letter [ 1 ], we studied the chi- 
ral properties of an SU(3) gauge theory with Nf massless 
Dirac fermions in the fundamental representation as Nf is 
increased from 2 to 6. We noted that the Nf = 2 simu- 
lations are in good agreement with measured QCD values, 
and that the Nf = 6 results indicate substantial enhance- 
ment of the chiral condensate. Here we extend our study of 
these two theories, presenting results for the electroweak S 
parameter and for the lightest vector and axial resonances. 

There is evidence from lattice simulations [|2}^6) that for 
an SU(iV) gauge theory, infrared conformality exists for 
a range of Nf values from the onset of asymptotic free- 
dom down to a critical value Nf. Below this "conformal 
window", chiral symmetry breaking and confinement set 
in. Even for Nf < Nf there can remain an approximate 
infrared fixed point provided that < Nf — Nf <C Nf. 
The scale of chiral symmetry breaking is then small, and 
the fixed point approximately governs the theory from the 
breaking scale out to some higher scale. This "walking" 
behavior leads to chiral-condensate enhancement, which 
can address the problem of obtaining large enough quark 
and lepton masses in technicolor theories. 

It has been suggested 0^9) that walking theories could 
address another problem by leading to smaller values of the 
electroweak S parameter. The value of S is related to the 
spectrum of vector and axial resonances in the theory. As 
in Ref. [T), we start with Nf = 2, allowing us to check the 
reliability of our methods by comparison with QCD phe- 
nomenology. Proceeding carefully toward Nf is prudent 
since the emergence of widely separated scales associated 
with walking is challenging for lattice methods. 

We first compute the S parameter from the defining cur- 
rent correlators, and then present results for the lowest ly- 
ing vector and axial masses and decay constants. We dis- 
cuss our results along with the related Weinberg spectral 
function sum rules, and then summarize. 



The S parameter The S parameter can be defined in 
terms of the vector and axial current-correlation functions 
with, by convention, the would-be Nambu-Goldstone- 
boson (NGB) contribution to the standard-model (SM) ra- 
diative corrections removed. With Nj/2 massless elec- 
troweak doublets, it can be written as ifTOl 
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where Ii vv (Q 2 ) and Haa(Q 2 ) are the transverse corre- 
lation functions for a single electroweak doublet, R(s) = 
12-7T Im n'(s), and m H is the reference Higgs mass. The 
presence of Ry(s) — Ra(s) in the spectral integral sug- 
gests that S could decrease if the resonance spectrum be- 
comes more parity doubled with increasing Nf. 

For Nf = 2, there are 3 Goldstone bosons, with the 
I 3 = 1 pair leading to Ry(s) — > 1/4 as s — > 0. 
(Ra(s) — > 0.) The standard-model subtraction removes 
the resultant infrared divergence. With Nf/2 electroweak 
doublets, there are Nf — 1 Goldstone bosons in the ab- 
sence of other interactions. Among these, (Nf/2) 2 pairs 
contribute to S, leading to R v (0) = (l/4)N f /2. With 
standard-model and other interactions included, Nf — 4 of 
the Nf — 1 Goldstone bosons will be pseudo-Goldstone 
bosons (PNGBs). The S parameter is again infrared finite, 
depending logarithmically on the masses of the PNGBs. 

Lattice simulations are carried out with a finite fermion 
mass rrif, requiring extrapolation to reach the chiral limit. 
With Nf/2 electroweak doublets, since we do not include 
SM and other interactions to give mass to the PNGBs, the 
extrapolation for Nf 7^ 2 would lead to log m / terms in S. 
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For our Nf = 6 simulations, m/ is not yet small enough 
to see clear evidence for these chiral logs. For smaller m/, 
the log m f terms would be replaced by logarithmic depen- 
dence on the PNGB masses in the full theory. 

Simulation Details Simulations are performed using 
domain-wall fermions and the Iwasaki improved gauge ac- 
tion iflTll . The domain- wall formulation suppresses the chi- 
ral symmetry breaking associated with fermion discretiza- 
tion, and preserves flavor symmetry at finite lattice spac- 
ing, both desirable properties for computation of the S- 
parameter. Gauge configurations are generated as in Ref. 
I2. Dimensionful quantities are given in lattice units. 

The lattice volume is set to 32 3 x 64, with the length of 
the fifth dimension L s = 16 and the domain-wall height 
m = 1.8. The choices (3 = 2.70 for N f = 2 and P = 
2.10 for Nf = 6 lead to nearly the same physical scale in 
lattice units. Simulations are performed for fermion masses 
m f = 0.005 to 0.03, although the N f = 2 results for 
m/ = 0.005 may suffer from finite- volume effects, and 
are not included in the analysis. At finite lattice spacing, 
even with m,f = 0, the chiral symmetry is not exact, with 
the violation captured in a residual mass m res <C mf. The 
total fermion mass m is then m = rrif + m res . 

Current Correlators The lattice expression for the cur- 
rent correlator of interest is 



^vv(Q) 



5^U VV (Q 2 ) - (Q»Q"/Q 2 )Il vv (Q 2 ) 
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and similarly for Ii AA . Here V M is the conserved domain- 
wall vector current, V v is the non-conserved local cur- 
rent, and Z is a non- perturb ative renormalization constant. 
(x + /t/2) appears because V M (x) is point split on the 
link (x,x + /_i). The use of conserved currents ensures 
that lattice artifacts cancel in the V — A current correlator 

n v - A {Q 2 ) = n vv (Q 2 ) - n AA (Q 2 ) E2. 

We calculate Uv-a(Q 2 ) f° r a range of positive (space- 
like) Q 2 values, and for each m t extrapolate to Q 2 = to 
determine the slope 4:irU' v _ A (0) entering the S parame- 
ter. In Fig. [I] we show the simulation data for U V ^ A (Q 2 ), 
along with fit curves. The data itself indicates that for 
Nf = 2, H' V _ A {0) increases at smaller rrt/ values, while 
for Nf = 6, it decreases, already suggesting a relative de- 
crease in S per electroweak doublet at Nf = 6. We fit 
the Hv~a(Q 2 ) data for Q 2 < 0.4 using a four-parameter, 
Pade(l,2) form (linear numerator, quadratic denominator). 
These fits, behaving like l/Q 2 at large positive Q 2 , are 
shown with statistical error bands in Fig. [T] Each has two 
poles at real, negative Q 2 , but they represent a time-like 
structure with cuts and multiple poles. Each fit leads to 
a value of H' V _ A (0) stable as the number of Q 2 points is 
varied. 

The correlator slopes at Q 2 = are plotted in Fig. 2. In 
this figure and others to follow, we plot versus Mp/My 
rather than m, where M P is the Goldstone-boson mass Q, 
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FIG. 1: Il v _ A (Q 2 ) data and fits for Nf = 2 and 6. Fits, over the 
range Q 2 < 0.40, are done separately for each m/. 
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FIG. 2: V - A correlator slopes at Q 2 = for Nf = 2 (red dia- 
monds) and Nf — 6 (blue circles). For each of the solid points, 

M P L > 4. 



and M V o is the extrapolated mass of the lightest vector 
state. We plot in this way since the relation between Mp 
and m is strongly iVy -dependent. The value of M vo , to 
be discussed later, is roughly 0.2 in lattice units for both 
Nf = 2 and 6. For each Nf = 6 point and for the five 
heaviest Nf = 2 points, MpL > 4, keeping the pion 
Compton wavelength well inside the lattice. 

As anticipated from inspection of the data in Fig. [TJ 
n' y _ j4 (0) at N f = 6 drops below W V _ A (0) at N f = 2 for 
the smaller Mp values, suggesting a suppression of S at 
Nf = 6. This interpretation requires care, however, since 
the extrapolation Mp oc m — > is dominated by chiral 
logs (~ log(l/m)) for both Nf = 2 and 6. 
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S-Parameter Results The S parameter (Eq. [TJ is sim- 
ply the correlator slope multiplied by the number of elec- 
troweak doublets, with the SM subtraction. We estimate 
the SM subtraction by evaluating the ASsm integral in 
Eq. [I] with an infrared cutoff at s = 4Mp, and taking 
rriH = M vo . For the case 2M P < M vo , 
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We use values for M P and M vo determined in Ref. Q- 
The choice m H = M vo corresponds roughly to a 1 TeV 
value for the reference Higgs mass. 

In Fig. [3) we plot S = 4vr(iV / /2)n' y _ A (0) - AS SM - 



For Nf = 2, the results are consistent with previous lattice 
simulations [El [13). The SM subtraction at N f = 2 is 
small, reaching a value ~ 0.04 for the lowest solid mass 
point, corresponding to my = 0.010. A smooth extrapo- 
lation to m = is expected since the LO chiral logs even- 
tually appearing in H' V _ A (0) are canceled by the SM sub- 
traction, Eq. [3] Given the linearity and small slope of the 
solid data points, we include a linear fit and extrapolation. 
An NLO term of the form MpXogMp has not been ruled 
out, but it is not visible in our data. The fit, with error band, 
is shown in Fig. [5] giving S m = = 0.35(6), consistent with 
the value obtained using scaled-up QCD data flOll . 



The Nf = 6 results for S are also shown in Fig. [3] The 
SM subtraction is again very small as at Nf = 2. The 
important feature is that the value of S at the lower mass 
points drops below a value obtained by simply multiplying 
the Nf = 2 result by a factor of 3. (For an Nf = 6 theory 
with only a single electroweak doublet, the value of S at 
the lower rrif values of Fig. [3] would be well below that 
of the Nf = 2 theory.) This trend has set in at Nt = 6 
even though 6 <C Nf. As m is decreased further at Nf = 
6, S as computed here will eventually turn up since the 
SM subtraction leaves a chiral-log contribution. For Nf/2 
electroweak doublets, S ~ (1 / 12tt)[N] / A - 1] logMjL 
In a realistic context, the PNGBs receive mass from SM 



FIG. 4: Axial and vector masses, Ma and My, and their ratio. 
Straight lines show linear fits to the solid points (MpL > 4), with 
the extrapolated values and errors shown to the left. 



and other interactions not included here, and these masses 
provide the infrared cutoff in the logs. 

Vector and Axial Masses A question of general inter- 
est for an SU(N) gauge theory is the form of the reso- 
nance spectrum as Nf is increased toward Nf. A trend to- 
ward parity doubling, for example, would provide a strik- 
ing contrast with a QCD-like theory. If the gauge theory 
plays a role in electroweak symmetry breaking, then this 
trend could be associated with a diminished S parameter. 

We have so far computed the masses, M v and Ma, and 
decay constants, F v and Fa, of the lowest-lying vector and 
axial resonances. We plot the masses along with their ratio 
in Fig. [4j Since the data points for each case except Ma 
at Nf = 6 are quite linear, with a small slope, and since 
in each case, the NLO term in chiral perturbation theory 
is linear in Mp oc m, we include a linear fit to the solid 
points (MpL > 4). The error bars on the extrapolations 
are also shown. M v extrapolates to 0.215(3) for Nf = 2, 
and to 0.209(3) for N f = 6. 

For Nf = 2, the extrapolated value of Ma/M v = 
1.476(40) is roughly consistent with the experimental re- 
sult of 1.585(52) fll). The N f = 6 data points for M A 
do not yet allow a simple fit and extrapolation, but they 
do indicate a substantial decrease in Ma/Mv in the chiral 
limit. This trend toward parity doubling suggests that the 
spectrum could become even more parity doubled as Nf is 
increased further, toward Nj. 

Vector and Axial Decay Constants Our simulation re- 
sults for F v and F A are shown in Fig. [3] using the nor- 
malization conventions of Ref. ifTOl . The dependence on 
Mp /M vo is mild, and once again, for each case except the 
A at Nf = 6, quite linear with a small slope. Although 
there is known to be an NLO chiral log for the decay con- 
stants, it is not visible in the linear points, so we include 
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FIG. 5: Decay constants Fy and Fa. The straight lines represent 
linear fits to the Nf — 2 data, with the extrapolated values in the 
chiral limit shown at the left end of the lines. 



a linear fit to the Nf = 2 data. The linearly extrapolated 
values, converted to physical units using the lattice scale 
determined from M vo , are F v = 141.8(3.8) MeV and 
Fa = 138.9(8.2) MeV, agreeing well with the measured 
OCD results ifTTlfBl . 

Discussion The relation between a diminished S pa- 
rameter and the spectrum can be explored through the dis- 
persion relation 

Q 2 f°° ds R v (s) - R A (s) 
12tt 



n v -A(Q 2 



7T 
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s + Q 2 

where F P ~ 93 MeV for two-flavor QCD. Here, we in- 
clude a few remarks about this use of the dispersion rela- 
tion via the assumption of single-pole dominance 

As a preliminary, recall that in the chiral limit (and 
in the continuum), the operator product expansion gives 
U V -a(Q 2 ) ~ 1/Q 4 as Q 2 — > oo. In the resultant first 
and second Weinberg sum rules (WSRs), the assumption 
of single pole dominance, R v _ A {s) = Yin 2 F VA 5(s — 
M ViA ), then leads to the relations F v - F\ = Fj\ and 
F 2 M 2 - F\M\ = 0. At finite m, however, ILy_ A (Q 2 ) 
falls more slowly with Q 2 . The simulation data in the Q 2 
range of Fig. [T]behaves like m(ip'ip) / Q 2 , but there will be 
an additional small (0(m 2 )) term [16|. While the data are 
consistent with the first WSR, the single-pole-dominance 
relation Fy — F\ = Fp is not satisfied at Nf = 2 for any 
m, the LHS being less than 50% of the RHS. The extrapo- 
lated values and the QCD experimental values fail by even 
more. A similar discussion applies to the second WSR at 
Nf = 2 and the resultant pole-dominance relation. For 
Nf = 6, the errors on Fa do not yet allow a useful test. 

A crude expression for S, which derives from the Q 2 —> 
limit of the dispersion relation, can be obtained by assum- 
ing single-pole dominance, and neglecting the SM subtrac- 
tion. The result, S = 4tt[F v /M v - F\/M\], pays no 
heed to a reference Higgs mass and has no direct depen- 
dence on the PNGB masses for Nf = 6. On the other 
hand, the stronger UV convergence of the integral expres- 
sion (Eq. [TJ could make single -pole dominance more re- 



liable than for the WSRs. In fact, evaluation of the above 
expression using our simulation data leads, for Nf = 2, 
to a result within 30% of our direct simulation of S at the 
smaller Mp values. For Nf = 6, the agreement with the 
direct simulation is at least as good, and also shows the 
relative decrease of S per electroweak doublet. 

Conclusions We have described a lattice simulation of 
the masses and decay constants of the lowest-lying vec- 
tor and axial resonances, and the electroweak S parameter, 
in an SU (3) gauge theory with Nf = 2 and 6 fermions 
in the fundamental representation. The spectrum becomes 
more parity doubled and the S parameter per electroweak 
doublet decreases when Nf is increased from 2 to 6. The 
final value of S for any Nf > 2 will depend logarithmi- 
cally on the masses of PNGBs, generated by SM and other 
physics not included here. It will be interesting to study 
these trends as Nf is increased further, toward N*. 

Acknowledgements We thank the Aspen Center for 
Physics, and the LLNL Multiprogrammatic and Institu- 
tional Computing program for time on the BlueGene/L su- 
percomputer. This work was supported by the NNSA and 
Office of Science of the U.S. Department of Energy, and 
by the U.S. National Science Foundation. 



[1] T. Appelquist, A. Avakian, R. Babich, R. C. Brower, 
M. Cheng, M. A. Clark, S. D. Cohen, G. T. Fleming, 
J. Kiskis, E. T. Neil, J. C. Osborn, C. Rebbi, D. Schaich, and 
P. Vranas, Phys. Rev. Lett. 104, 071601 (2010), 0910.2224. 

[2] T. Appelquist, G. T. Fleming, and E. T. Neil, Phys. Rev. Lett. 
100, 171607 (2008), 0712.0609; Phys. Rev. D79, 076010 
(2009), 0901.3766. 

[3] Z. Fodor, K. Holland, J. Kuti, D. Nogradi, and C. Schroeder 
(2009), 0907.4562. 

[4] K.-i. Nagai, G. Carrillo-Ruiz, G. Koleva, and R. Lewis 
(2009), 0908.0166. 

[5] L. Del Debbio, B. Lucini, A. Patella, C. Pica, and A. Rago, 
Phys. Rev. D82, 014510 (2010), 1004.3206. 

[6] A. Hasenfratz, Phys. Rev. D82, 014506 (2010), 1004.1004. 

[7] T. Appelquist and F. Sannino, Phys. Rev. D59, 067702 
(1999), hep-ph/9806409. 

[8] S. D. H. Hsu, F. Sannino, and J. Schechter, Phys. Lett. B427, 
300 (1998), hep-th/9801097. 

[9] M. Kurachi and R. Shrock, Phys. Rev. D74, 056003 (2006), 
hep-ph/0607231. 
[10] M. E. Peskin and T. Takeuchi, Phys. Rev. D46, 381 (1992). 
[11] C. Allton et al. (RBC-UKQCD), Phys. Rev. D78, 114509 

(2008), 0804.0473. 
[12] P. Boyle, L. Del Debbio, J. Wennekers, and J. Zanotti (RBC- 
UKQCD), Phys. Rev. D81, 014504 (2010), 0909.4931. 
[13] E. Shintani et al. (JLQCD), Phys. Rev. Lett. 101, 242001 

(2008), 0806.4222. 
[14] C. Amsler et al. (PDG), Phys. Lett. B667, 1 (2008). 
[15] N. Isgur, C. Morningstar, and C. Reader, Phys. Rev. D39, 
1357 (1989). 

[16] M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Nucl. 
Phys. B147, 448 (1979). 



